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^P , Abstract: Dispersive representations of the Kir vector and scalar form factors are used to 

fit the spectrum of r — )■ Kiruj- obtained by the Belle collaboration incorporating constraints 
from results for Ki^ decays. The slope and curvature of the vector form factor are obtained 

y^ ^ directly from the data through the use of a three-times-subtracted dispersion relation. 

sS ■ We find A'+ = (25.49 ± 0.31) x 10"^ and \'[ = (12.22 ± 0.14) x 10"^. From the pole 

position on the second Riemann sheet the mass and width of the i^*(892) are found 
to be w._ft-*(892)± = 892.0 ± 0.5 MeV and '^K*{m2)^ = 46.5 it 1.1 MeV. The phase-space 
integrals needed for i^/3 decays are calculated as well. Furthermore, the Kt: isospin-1/2 
P-wave threshold parameters are derived from the phase of the vector form factor. For the 
scattering length and the effective range we find respectively a^ = (0.166 it 0.004) m~^ 
and h]/"^ = (0.258 ± 0.009) m'^ 
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1. Introduction 

The differential decay distributions of X — )• vr / f/ {Ki.^ ) and r — )■ Kiri^r decays are governed 
by two Lorentz-invariant Kir form factors that encode the non- perturb at ive physics, namely 
the vector, denoted F^'^{q'^), and the scalar, F^'"{q'^). According to the kinematical con- 
figuration, q represents the exchanged (Ki^) or the total (r — )• Ktti/t-) Ktt four- momentum. 
A good knowledge of these form factors is of fundamental importance for the determination 
of many parameters of the Standard Model, such as the quark-mixing matrix element \Vus\ 
obtained from Ki^ decays [1], or the strange-quark mass mg determined from the scalar 
QCD strange spectral function [2]. Recently, several collaborations have produced data 
for Ki^ decays and new high-statistics data for r — )• KttVt have been published by the B 
factories. The new data sets provide the substrate for up-to-date theoretical analyses of 
the Ktt form factors. 

Historically, the main source of experimental information on Ktt form factors have been 
Ki^ decays. Recently, five experiments have collected data on semileptonic and leptonic 
K decays: BNL-E865 [3], KLOE [4], KTeV [5], ISTRA+ [6], and NA48 [7]. The resuhs 
from these analyses yielded an important amount of information on form factors as well as 
stringent tests of QCD at low-energies and of the Standard Model itself (for recent reviews 
on theoretical and experimental aspects of kaon physics we refer to Refs. [8,9]). Additional 
knowledge on the Ktt form factors can be gained from the dominant Cabibbo-suppressed 



r decay: the channel r — >• KirVr- The r is the only known lepton heavy enough to decay 
into hadrons and its hadronic decays constitute a rather clean environment for the study of 
QCD at relatively low energies [10] and notably for the determination of the QCD coupling 
Us [11-14]. In the 1990s, the i^vr spectrum for r — )• KnVr was measured by ALEPH [15] 
and OPAL [16]. Lately, however, the B factories have become a superior source of high- 
statistics data for this reaction by virtue of the important cross-section for e"^e~ — )• t^t^ 
around the T(4S') peak. As a result, as many as 10^ r pairs were recorded by Belle and 
BaBar [17]. A detailed spectrum for r — )• Ksn^Ur produced and analysed by Belle was 
published in 2008 [18] with an event sample larger than in the LEP experiments by almost 
a factor of 65, allowing for a detailed analysis of its shape. Also, a preliminary BaBar 
spectrum with similar statistics has appeared recently in conference proceedings [19] and, 
finally, BESIII should produce results for this decay in the future [20]. 

On the theory side, a salient feature of the form factors in the kinematical region 
relevant for Ki^ decays, i.e. mf < q^ < {mx — fnirY'-, is that they are real. Within 
the allowed phase-space they admit a Taylor expansion and the energy dependence is 
customarily translated into constants A^q defined as^ 



F+fl{q^) = F+,o(0) 






fl.l) 



In r — )■ KttVt decays, however, since {niK + w-tt)^ < q^ < m^, one deals with a different 
kinematical regime in which the form factors develop imaginary parts, rendering the ex- 
pansion of Eq. (|1 . 1| ) inadmissible. One must then resort to more sophisticated treatments. 
Moreover, in order to fully benefit from the available experimental data, it is desirable 
to employ representations of the form factors that are valid for both Ki^ and r — )• KttVt 
decays. In Ref. [21], a new expression for i*+(s) was derived within the Resonance Chiral 
Theory (RChT) framework [22] and, subsequently, the authors reanalysed the Belle spec- 
trum for the decay r — )• Kiriyr with success [23]. This analysis yielded new values for the 
constants A'_^ and A'_^ emphasising the interplay between r — )• KttVt and Ki^ experiments. 
Prom general principles of analyticity, the form factors must fulfil a dispersion relation. 
Unitarity provides an additional constraint on the imaginary part of the form factors, 
rendering possible the design of dispersive representations of F+ and Fq that are suited to 
describe both r — )■ KnUr and Ki^ decays. For the vector form factor, a step towards this feat 
was taken in Refs. [24,25] where we introduced several subtracted dispersive representations 
of -F+. Our final proposal was a three-times-subtracted dispersive representation in which 
A'_|_ and X'\_ are parameters that were determined via a successful fit to the Belle spectrum. 
A similar dispersive approach to F-f was presented in Ref. [26] and has been used by the 
KTeV collaboration to fit their Ki,^ spectra [27]. Finally, a dispersive representation for 
-F+ that includes inelastic effects was introduced in Ref. [28]. Concerning the scalar form 
factor, a thorough description that takes into account analyticity, unitarity, the large- A'^c 
limit of QCD, and the coupling to Krj and Rrj' channels was introduced in Ref. [29] and 
updated in Refs. [2,30,31]. Another single-channel dispersive representation of Fq can be 



""^From now on we refrain from writing the superscript Kiv on tiie form factors. 



found in Ref. [32] and a description based on the so-called method of unitarity bounds was 
recently presented in Ref. [33]. 

The main purpose of our paper is to produce an analysis of the Belle spectrum for 
r — )• KttVt incorporating constraints from experimental results on Ki^ decays. We have 
already advocated that such a combined treatment of both reactions could further our 
knowledge of the form factors hence paving the way for a better determination of \Vus\ [25]. 
Moreover, we aim at extracting as much information as possible from the r — t- KttUj- 
spectrum. With the present statistics the spectrum allows for a study of Ktt dynamics in 
the P wave, which gives the prevailing contribution to the decay. Watson's theorem [34] 
guarantees that below inelastic thresholds the phase of the form factor equals the scattering 
phase and, therefore, one can perform a study of the dominant Ktt P-wave threshold 
parameters. In addition, it has been shown [23, 24] that the present statistics permits a 
competitive determination of the pole position of the K*(892) as well as the position of 
a second vector resonance, although less precisely in the latter case. Here, we determine 
these two poles exploiting a novel strategy in which fits are done directly in terms of the 
physical pole positions on the second Riemann sheet. This improvement with respect to 
previous works [18,23,24] yields a determination of the pole positions with a better control 
of uncertainties and correlations. 

In our analysis, for the vector form factor we employ the dispersive representation of 
Ref. [24] whereas for the scalar Ktt form factor we use the up-to-date results of Ref. [31]. 
Since the details of these descriptions can be found in the original works, here we shall 
concentrate on the results that arise from our fit, namely i) the pole positions for the 
K*{892)^ and i^*(1410)^ resonances, ii) A^ and A" , in) the result of the phase-space 
integrals needed in Ki^ decays, and iv) the Kit isospin-1/2 P-wave scattering phase and 
the respective threshold parameters. 

Our paper is organised as follows. First, in Sec. |^, we briefly review the dispersive 
treatment of the vector and scalar Kir form factors. Then, in Sec. |3|, we present a fit to 
r — )• KiTVr data alone. In Sec. H, the results for a fit incorporating constraints from Ki^ 
experiments are given and, in Sec. 13, we derive our results for the phase-space integrals 
relevant for Ki^ experiments. We discuss the results for the Ktt threshold parameters and 
scattering phase shifts in Sec. |6|. Our final results and a comparison with other results 
found in the literature are presented in Sec. 0. 

2. Dispersive Kiz form factors 

The Ktt form factors are defined as follows [8] 



1 1 

m\r — mt 



{k + vY-^^hr^ik-py 



Q J Q 

(2.1) 

where F^{q'^) and Fo(q^) are the vector and scalar form factors respectively and q^ = 

(k — p)"^. It follows from the definition that both form factors share the same normalisation 

at zero i*+(0) = Fq{0). For convenience, we work with normalised form factors F-^-,o{q'^) 



such that 

#«tf).^. (2.2, 

First, in determinations of \Vus\, a rehable value for the normahsation at zero is crucial in 
order to disentangle the product |V^s|F+(0). In this respect, Chiral Perturbation Theory 
(ChPT) and lattice QCD are the most trustworthy methods to obtain F+(0). Here we 
are concerned with another aspect of the form factors, namely their energy dependence 
encoded in F^fi{q'^). The precise knowledge of F^fi{q'^) is needed when performing the 
phase space integrals for Ki^ decays or when studying in detail the r — )• K-KUr spectrum. 
Finally, one should bear in mind that when considering r decays, one deals with a crossing- 



2^ 



symmetric version of Eq. (2.1) for the Ktt pair is in the final state. In this case, q = s 
{k + p)"^ > {rriK + ttItt)'^ and the form factors develop imaginary parts. 

In T — 7- K-kUt and K^,^ decays, the term containing the vector form factor F^(q 
dominates the differential decay widths. The form factor, in its turn, receives a prevailing 
contribution from the K*{892). This fact motivated the description of Refs. [21,23] within 
RChT, which was based on an analogous treatment of the pion vector form factor [37,38]. 
Although dominated by the K*{892), the authors of Refs. [21,23] noted that a second 
resonance, identified with the i^*(1410), must be included in F^{s) to account for the 
higher-energy part of the r — )• KTrVr spectrum. The description of Refs. [21,23], albeit 
successful, has a slight drawback, namely it satisfies the analyticity constraints only in a 
perturbative sense. Although the violation of analyticity is expected to be of higher orders 
in the chiral expansion, a description based on a dispersive treatment was necessary to 
corroborate this pattern. In Ref. [24] we designed such dispersive representations of -F+(s). 

The rationale for our approach is as follows. From general principles, the form fac- 
tor must satisfy a dispersion relation. Supplementing this constraint with unitarity, the 
dispersion relation has a well-known closed-form solution within the elastic approximation 
referred to as the Omnes representation [39]. Although simple, this solution requires the 
detailed knowledge of the phase of i*+(s) up to infinity, which is unrealistic. An advanta- 
geous strategy to circumvent this problem is the use of additional subtractions, as done for 
the pion form factor in Ref. [40]. Subtractions in the dispersion relation entail a suppression 
of the integrand in the dispersion integral for higher energies. An n-times-subtracted form 
factor exhibits a suppression of s^'^"'^^' in the integrand. Thereby, the information that 
was previously contained in the high-energy part of the integral is translated into n — 1 
subtraction constants. In Ref. [24] we performed fits to the Belle spectrum of r — )• KnUr 
varying the number of subtractions and testing the description with one and two vector 
resonances. The outcome of these tests, described in detail in Ref. [24], is that for our 
purposes an optimal description of F^{s) was reached with three subtractions and two 
resonances. Here we quote the resulting expression 



F+{s) = exp 



ai 



+ 7:«2^ — \ ds' 



(2.3) 



m^_ 2 m^_ vr J {s')^{s' — s — iO) 

In the last equation, sk-k = ("^i^o + ?tt.^-)^ and the two subtraction constants ai and a2 are 
related to the Taylor expansion of Eq. (|1.1D as A^ = qi and A'_^ = a2 + Oi\. It is opportune 



to treat them as free parameters that capture our ignorance of the higher energy part 
of the integral. The constants A^ and A^ can then be determined through the fit. The 
main advantage of this procedure, advocated for example in Refs. [24,26,32,40], is that 
the subtraction constants turn out to be less model dependent as they are determined by 
the best fit to the data. The calculation of these constants, on the other hand, depends 
strongly on the perfect knowledge of S{s). However, since now ai^2 are determined by the 
data, in the limit s — )• oo the asymptotic behaviour of -F+(s) cannot be satisfied. This is so 
because a perfect cancellation between terms containing qi and 02 with polynomial terms 
coming from the dispersion integral must occur in order to guarantee that -F+(s) vanishes 
as 1/s. We have checked that our form factor, within the entire range where we apply it 
(and beyond), is indeed a decreasing function of s which renders this approach credible. 

With Eq. ( |2.3| ) , the transition from the kinematical region of r — t- KirUr to that of Ki^ 
decays is straightforward and the dominant low-energy behaviour of -F+(s) is encoded in A^ 
and A'_|_ . The cut-off Scut in the dispersion integral is introduced to quantify the suppression 
of the higher energy part of the integrand. The stability of the results is checked varying this 
cut-off in a wide range from 1.8 GeV < ^Scut < 00. It is important to stress that Eq. ( |2.3| ) 
remains valid beyond the elastic approximation provided 6{s) is the phase of the form 
factor, instead of the corresponding scattering phase. But, of course, in order to employ it 
in practice we must have a model for the phase. As described in detail in Ref. [24] , we take a 
form inspired by the RChT treatment of Refs. [21,23] with two vector resonances. Here we 
relegate the details concerning 6 to Appendix^. However, one important remark is in order. 
Since we keep the real part of the loop bubble integral Re^(s) in Eq. (|A.3| ), the mass and 
width parameters of Ref. [24] are shifted as compared with those of Refs. [18,21,23]. This 
shift emphasises the need for the computation of the physical pole position of the resonances 
on the second Riemann sheet. We have shown [24, 25] that although the mass and width 
parameters from Refs. [18, 23, 24] differ considerably, the pole positions arising from the 
models are in good agreement. To clarify this issue further, in this work we implement a 
numerical improvement in our codes that allows us to perform the fits directly in terms of 
the pole positions on the second Riemann sheet. This new procedure is clearer as it avoids 
the cumbersome intermediate stage where one must compute the pole positions from the 
unphysical parameters to obtain meaningful results [36]. Furthermore, correlations and 
uncertainties are obtained directly for the physical poles and are therefore more reliable. 

In the previous analysis of Refs. [23, 24] the scalar form factor was shown to play an 
important role for the low-energy part of the r — )• KttUt spectrum, between threshold and 
~ 0.8 GeV. On the other hand, the fit was not very sensitive to the details of Fq as it is in 
the case of -F+. Therefore, we again rely on the coupled channel representation of Fq first 
presented in Ref. [29] and updated in Refs. [2,30,31]. The main features of this treatment 
can be found in Appendix ^ 

3. Fit to T -)■ Ktzv^ 

Before proceeding to a fit that combines information from r — t- KttUt and Ki^ data, we 
shall perform in this section a short update of Ref. [24]. The aim is twofold. First we want 



to ascertain the impact of performing the fit directly in terms of the physical pole positions 
for the vector resonances. Second, the results of this section serve as a point of reference 
for the new analysis. For the sake of completeness, we recall here how the Kit form factors 
enter the description of r — ;■ KnUr ■ 

Assuming isospin invariance, the differential decay distribution for r — t- Knuj- can be 
cast in terms of the Kir form factors as 
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where we summed over the two possible decay channels r — t- K^tt v^ and r -^ K ix^v^ 



that contribute in the ratio 2:1. In Eq. (3.1), S'ew is an electroweak correction factor. 



and 
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(k+p) with k and p being respectively the kaon and pion momenta. 
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(3.2) 



is the kaon momentum in the rest frame of the hadronic system. In order to analyse the 
data, one must rely on an ansatz for the number of events observed in a given bin of the 
experimental spectrum. As explained in Ref. [23] the theoretical number of events A^* in 
the i-th bin is taken to be 



N,. 



th 



■Nt\-^ 



dr^TT 



2 3 ""TrBK^ d^i 



(4), 



(3.3) 



where Mt is the total number of events, the factor ^ and 3 account for the fact that the 
KsT^~ channel was analysed, A^ is the width of the i-th bin, F^ is the total r decay width, 
Bktt is a normalisation constant that, for a perfect description of the spectrum, should be 
the T —7- Kt:Ut- branching ratio, and, finally, s\^ is the centre of the z-th bin. In the case of 
Belle's spectrum [18] one has Nt = 53110 and a constant bin width Ab = 11.5 MeV. 
In this fit, we minimise the \^ function given by 
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(3.4) 



where N^^^ and a^exp are, respectively, the experimental number of events and the cor- 
responding uncertainty in the i-th bin. The prime in the symbol of sum indicates that 
bins 5, 6, and 7 are excluded from the minimisation^. In the x^i following a suggestion 
of the experimentalists [35], we include data up to bin number 90 which corresponds to 



^If these three points are included in the fit the results do not change significantly although the x^ is 
larger. Furthermore, there is no indication for a peak at this energy and BaBar spectra do not display a 
bump close to threshold. Hence, we decided, following Refs. [23,24], to exclude these points. For a visual 
account, points not included in the x^ are shown as unfilled circles in Fig. hi 



y/s = 1.65925 GeV. Finally, the lowest data point is not taken into account since, with 
physical meson masses, its centre lies below the Ktt threshold. The second term on the 
right-hand side of Eq. ( |3.4| ) was not included in the x^ function of Ref. [24]. It introduces 
an additional restriction that allows us to treat the normalisation -B/^tt of Eq. ( |3.3D as a 
free parameter. Then, the parameters of the fit are 8 in total. First, the two constants A^ 
and A" responsible for the behaviour of i^+(s) near the origin. Second, the five parameters 
that determine the resonance properties, i.e. the complex pole positions of the i^*(892) 
and'^ i^*(1410) and the mixing parameter 7 [see Eq. (|A.2D ]. The mass and width of the 
resonances are extracted from the complex pole position sr as [36] 

1 

^/s^ = mR- -Tr. (3.5) 

The phase of the form factor is fully determined by the latter set of parameters. The 8th 
parameter of the fit is the normalisation Bxn- 

In the fit, we employ the following numerical values: |T4s|-^+(0) = 0.2163(5) [9], Gp = 
1.16637(1) X 10-^ GeV-2 [41], m^ = 1776.84 MeV [41], Sew = 1.0201(3) [42], /^ = 
92.21(14) MeV [43], fx/U = 1-197(6) [43], and 5^^ = 0.418(11)% [44,45]. We recaU that 
the cut-off Scut of Eq. (|2.3| ) has to be varied in order to check the stability of the results 
upon the high-energy part of the dispersion integral. When quoting final results one must 
therefore include an uncertainty due to the small residual dependence on Scut- The results 
of fits with four values of Scut, namely Scut = 3.24 GeV^, 4 GeV^, 9 GeV^, and Scut — ^ co, 
are displayed in Tab. ||. 





s^^t = 3.24 GeV^ 


Scut = 4 GeV2 


Scut = 9 GeV2 


Scut -^ 00 


BKn 


0.416 ± 0.011% 


0.417 ±0.011% 


0.418 ±0.011% 


0.418 ±0.011% 


(B't) 


(0.414%) 


(0.414%) 


(0.415%) 


(0.415%) 


niK* [MeV] 


892.00 ±0.19 


892.02 ±0.19 


892.03 ±0.19 


892.03 ±0.19 


Vk* [MeV] 


46.14 ±0.44 


46.20 ± 0.43 


46.25 ±0.42 


46.25 ± 0.42 


niK" [MeV] 


128111 


1280+1 


12781^^ 


12781^^ 


Tk*' [MeV] 


2431?^ 


19311 


1771^^ 


177l| 


7 X 10^ 


5 1+1-7 


-3.9li| 


_3 4+i-i 


4+1.1 


A'+ X 103 


24.15 ±0.72 


24.55 ± 0.68 


24.86 ±0.66 


24.88 ± 0.66 


a; X 10^ 


11.99 ±0.19 


11.95 ±0.19 


11.93 ±0.19 


11.93 ±0.19 


xVn.d.f. 


74.1/79 


75.7/79 


77.2/79 


77.3/79 



Table 1: Results for the fit to Belle's r — ?► Kirvr spectrum [18]. As a consistency check, for each 
one of the fits we give the value B^^^ obtained from the integration of Eq. (3.1). 



Some of the results of Tab. |l| are to be compared with those of Tab. 4.2 of Ref. [24]. 
Concerning A^, A" and 7 they are very similar if not identical. However, in Ref. [24], the x^ 
that was minimised did not include the second term in the right-hand side of Eq. ( |3.4D and 
therefore Bk-k was kept fixed or, otherwise, the strong positive correlation between Bktt 
and the constant A', would render a good determination of these parameters impracticable. 



'For simplicity, in tables we refer to the ii'*(892) and the A'* (1410) simply as K* and K" respectively. 



The main difference between tlie two fits lies, as already stressed, in the pole positions of the 
vector resonances. In Tab. |l[ the results correspond to physical masses and widths obtained 
from the pole positions in the second Riemann sheet according to Eq. ( |3.5| ) whereas the 
results of Tab. 4.2 of Ref. [24] are non-physical parameters. Consequently, results for 
masses and widths presented here should not be directly compared with the parameters 
of Tab. 4.2 of Ref. [24]. Instead, one should compare with the physical poles that can be 
found in Eqs. (5.4) and (5.5) of Ref. [24]. One sees that the central values agree nicely. On 
the other hand, the treatment of the uncertainties affecting the poles is here much more 
trustworthy. The results of Tab. come from an analysis performed by the MINOS function 
of the CERN-Minuit library. The errors are smaller than the ones quoted in Ref. [24] due 
to the proper inclusion of correlations. Finally, the fit is very stable against changes in 
Scut- This is specially true for the mass and width of the K*{892) but in all other cases 
variations are at most at the level of one standard deviation. In order to produce a feeling 
for the correlation coefficients between the parameters of our fits, as an example we display 
in Tab. |2| those corresponding to Scut = 4 GeV^. 





Bk. 


rriK* 


Tk* 


rriK*' 


Tk" 
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a; 


a: 


rriK* 


-0.119 


1 














Tk' 


0.041 


-0.017 


1 












niK*' 


-0.048 


-0.168 


-0.158 


1 










Tk*' 


0.110 


0.182 


0.303 


-0.628 


1 








1 


-0.148 


-0.244 


-0.425 


0.558 


-0.865 


1 






a; 


0.711 


0.008 


0.543 


-0.298 


0.462 


-0.653 


1 




a; 


0.880 


-0.132 


0.421 


-0.212 


0.355 


-0.466 


0.934 


1 



Table 2: Correlation coefficients for the parameters of the fit with Scut — 4 GeV^, third column of 
Tab. |. 



4. Fit to T — >■ KtvUt with restrictions from Ki^ 



Dispersive representations of Kit form factors can be used in order to simultaneously fit 
both r — )• Ktti't and Ki^ spectra. We have recently advocated [25], performing Monte Carlo 
(MC) simulations, that the main benefit of such a combined fit would be the reduction of 
the uncertainties on the parameters A^ and A^, leading to smaller uncertainties in the 
phase-space integrals needed for the extraction of Vus from kaon decays. For the want of 
an unfolded data set from Ki^ experiments, we perform here a fit to r — )■ KttVt- constrained 
by results for A^ and A" obtained from a compilation of Ki^ analyses [9]. 

In results obtained from quadratic representations such as the one of Eq. ( |l . l]) , the 

(n) 

errors on A^ q have a clear statistical meaning. In principle, therefore, it is straightforward 
to include that information in the x^ that is to be minimised by the fit. In this case, 
the statistical correlation between A^ and A^ must be taken into account. The x^ to be 



minimised contains then one additional term 






where the first two terms in the right-hand side are the same as in Eq. (3^) whereas the 
last one encodes the information from Ki^ analyses. In this last term, the vectors A^.''^^^ 
are given by 



. 'th,exp 

X': 



^+ — I ,//th,exp ) (4-2) 



+ / 

and the 2x2 matrix V is the experimental error matrix for A+ such that 

Vij = pij Gi aj , (4.3) 

where the indices refer to A^ and A" , pij is the correlation coefficient {pij = 1 if i = j), and 
(Tj the experimental errors on A^ and A" . For the experimental values we employ the results 
of the compilation of Kl analyses performed by Antonelli et al. for the FlaviaNet Working 
Group on Kaon Decays in Ref. [9]: X':^""^ = (24.9 ± 1.1) x lO^^ X'^^'^'P = (16 ± 5) x 10""^ 
and py ,A" = —0.95. Results for fits using the x^ function of Eq. ( p] ) with Scut = 3.24 
GeV^, 4 GeV^, 9 GeV^, and Scut — >■ oo are shown in Tab. |3|. In Fig. 0, the Behe spectrum 
for r — )• K-kVt is confronted with the results for the fit with Scut = 4 GeV^. Finally, as an 
example, the correlation matrix for Scut = 4 GeV^ is given in Tab. 0. 





Scut = 3.24 GeV^ 


Scut = 4 GeV2 


Scut = 9 GeV2 


Scut -^ OO 


Bk. 


0.429 ± 0.009 


0.427 ± 0.008% 


0.426 ± 0.008% 


0.426 ± 0.008% 


{B'D 


(0.426%) 


(0.425%) 


(0.423%) 


(0.423%) 


rriK* [MeV] 


892.04 ± 0.20 


892.02 ±0.20 


892.03 ±0.19 


892.03 ±0.19 


Tk* [MeV] 


46.58 ± 0.38 


46.52 ±0.38 


46.48 ± 0.38 


46.48 ± 0.38 


m^./ [MeV] 


1257130 


126811 


1270+24 


12711^4 


Tk*' [MeV] 


02-I+95 
"3^-^-76 


2381^^ 


2061^0^ 


2051^7 


7 X 10^ 


S 9+2.2 


r. 4+1.4 


-4.41}:^ 


_4 4+1-2 


A'+ X 103 


25.43 ± 0.30 


25.49 ±0.30 


25.55 ±0.30 


25.55 ±0.30 


a; X 10^ 


12.31 ±0.10 


12.20 ±0.10 


12.12 ±0.10 


12.12 ±0.10 


xVn.d.f. 


77.9/81 


78.1 /81 


79.0 /81 


79.1/81 



Table 3: Results of fits to Belle spectrum [18] of t — ^ KttVt with constraints from the Ki^ analysis 
of Ref. [9]. The x^ function is defined in Eq. (4.1). As a consistency check, for each one of the fits 



we give the value B^^ obtained from the integration of Eq. ([3.1|). 

Comparing the results of the fit constrained by Ki^ analyses. Tab. H, with the results 
of the fit to T — )■ KTrVr alone. Tab. 0, one sees that the statistical uncertainty in A^ and 
A^ is reduced roughly by a factor of 2. Another advantage of the new fit is that the 
results for A^ and A" are much more stable against changes in Scut- The errors in A^ are 
largely dominated by statistics in sharp contrast with Tab. |l] where the model dependent 
uncertainties arising from the Scut dependence were of the same order as the statistical 
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Figure 1: Fit result for the spectrum of r -^- Kirv^- witlr s^ut = 4 GeV^, third column of Tab. 0. 
The data are from the Belle collaboration [18]. Points represented with unfilled circles are excluded 
from the fit (see text in Section ||) . The solid red line represents the full fit including contributions 
from F+(s) and -fo(s)- The scalar contribution alone is represented by the dot-dashed orange line 
whereas the dashed blue line gives the vector contribution. 

ones. The central results of A^ and A" exhibit a small shift because Ki^ experiments 
favour larger values. The mass of the -ftr*(892), in its turn, turns out to be almost the same 
as in the previous fit and is still very stable with respect to changes in Scut- The ir*(892) 
width is slightly larger than before but compatible within one sigma with the previous 
result. The parameters of the second resonance have still large uncertainties but remain 
compatible with the results of Tab. |l|. Finally the normalisation Br.,^ turns out larger than 
in the previous fit due to a positive correlation with A^ and A" but fully compatible with 
the experimental experimental average B^^^ = 0.418(11). 





BKn 


ruK* 


Tk- 


rriK*' 


Tk*' 


7 


a; 


a; 


niK* 


-0.193 


1 














Tk* 


-0.414 


-0.007 


1 












rriK*' 


0.223 


-0.233 


-0.043 


1 










^K" 


-0.261 


0.243 


0.130 


-0.675 


1 








7 


0.399 


-0.344 


-0.193 


0.630 


-0.886 


1 






a; 


0.316 


0.058 


0.290 


-0.186 


0.252 


-0.386 


1 




a; 


0.776 


-0.233 


0.045 


-0.001 


0.054 


0.006 


0.747 


1 



Table 4: Correlation coefficients for the parameters of the fit with Scut — 4 GeV^, third column of 
Tab. U. 
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5. Ki^ phase space integrals 

From the results of our fits shown in Tab. |^ one can calculate the phase-space integral 
needed in the computation of Ki^ decay widths. The phase-space integral is defined as"^ 

{mK-m.T,)'^ 2 






where mi is the mass of the lepton and 



(5.1) 



\{t) = 1 + t^/m\^ + rl-2rl-2rl t/m\ - 2 t/m 



(5.2) 



with Ttt = rn^/rrij^. In the phase-space integral for the decays with an electron in the final 
state, Iks ! the smallness of the electron mass makes the contribution of Fq immaterial. The 
scalar form factor gives nevertheless a non-negligible contribution for X^g decays. In phase- 
space integrals for decays of charged kaons, we have assumed that the normalised form 
factors -F+,0 are isospin invariant, which amounts to assuming that isospin breaking effects 
are solely contained in F+(0). Then, for the phase-space factors of charged- kaon integrals 
we employ the mass of the charged kaon and that of the neutral pion. Tab. g contains our 
results for the integrals. In order to take into account all errors and correlations, a MC 
sample of parameter values employing the results from Tabs. |3| and ^ was generated. The 
integrals were computed for each set of parameters in these samples. The errors quoted in 
Tab. 1^ are of a gaussian nature to a good approximation. 





s,^t = 3.24 GeV^ 


Scut = 4 GeV2 


Scut = 9 GeV2 


Scut -^ 00 


^K2, 


0.15463(17) 


0.15465(16) 


0.15468(16) 


0.15468(16) 


^n. 


0.10275(10) 


0.10276(10) 


0.10277(10) 


0.10277(10) 


^K+ 


0.15900(17) 


0.15902(16) 


0.15905(16) 


0.15905(16) 




0.10573(11) 


0.10575(10) 


0.10576(10) 


0.10576(10) 



Table 5: Results for the phase-space integrals defined in Eq. ( p.l[ ) obtained with parameters from 
the fits of Tab. ||. The uncertainties include the statistical errors and correlations from the fit. 



6. Ktt isospin-1/2 P-wave scattering phase 

The decay r — t- KttUt offers a good environment for the study of Kn dynamics. From 
the point of view of strong interactions, the Ktt pair in the final state is isolated. As a 
matter of fact, this decay is certainly a better laboratory for the study of the Ktt phase 
than the hadronic reactions used in the classical determinations of the Ktt phase shifts. 



* We employ tlie notation of Ref. [8] but the definition of tlie integral is identical to that of Ref. [1] . 
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Figure 2: Phase of the form factor F+(s) together with experimental results from LASS [46] and 
Estabrooks et al. [47]. The opening of the first inelastic channel, K*Tr, is indicated by the dashed 
vertical line. The gray band represents the extrema from the fits of Tab. 0. 

Watson's theorem states that below the first inelastic threshold the form factors and the 
respective partial- wave scattering amplitudes share the same phase [34] . In the case of the 
P wave, the first inelastic channel one can consider is the quasi-two-body K*7r which opens 
at ~ 1030 MeV [28]. Therefore, below this value, the phase of our vector form factor can 
be compared with the respective scattering results. In Fig. ^, we compare our phase with 
those from LASS [46] and Estabrooks et al. [47]. In this comparison, one should bear in 
mind that isospin breaking effects could play a small role since the hadronic experiments 
measured the neutral channel whereas we have the charged one. Nevertheless, from the 
inspection of Fig. |2|, it is clear that our results are compatible with the experimental 
determinations of the Ktt I = 1/2 P-wave scattering phase shift between 850 MeV and 
roughly 1 GeV, just before inelasticity sets in. From threshold up to 850 MeV our results 
seem to be systematically lower than those from hadronic reactions. It is interesting to 
remark that the same behaviour is also observed in the recent Roy-Steiner-type analysis of 
Ktt scattering performed by Biittiker, Descotes-Genon and Moussallam [48]. Their phase 
is also somewhat below the experimental data up to about 950 MeV. Finally, we remind 
that the low-energy results from Estabrooks et al. [47] have been shown to be inconsistent 
with a dispersive analysis of Ktt scattering [49] and, unfortunately, LASS results [46] do 
not span the energy region close to threshold. 

Within the elastic domain, the phase of our form factor F-^-{s) equals the scattering 
phase for the P wave with 1=1/2. From the expansion of the corresponding partial wave 
i- matrix near threshold we can obtain the threshold parameters. Following Ref. [48], they 
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Scut = 3.24 GeV2 Scut = 4 GeV^ s^ut = 9 GeV^ Scut ^ oo 



mfTopTTo 0.1658(13) 0.1656(13) 0.1655(13) 0.1655(13) 

m^_ bl^"^ X 10^ 0.2573(24) 0.2581(23) 0.2582(23) 0.2583(23) 

ml_ c^'^ X 10^ 0.8987(81) 0.9001(76) 0.9000(75) 0.9000(75) 



Table 6: Threshold parameters defined in Eq. ( |6.lD calculated with the results of our main fit 
given in Tab. 0. The uncertainties are solely statistical. 

are defined for isospin / and angular momentum / as 

_|Ret/(,) = lsin2<5/(<?) = 5'^[af ^h\ q^ ^ c\ q' ^ 0{i)\ , (6.1) 

where g(s) is given by Eq. (^.2|). It is simple to express our results in the form of Eq. 
using 



s = m\ + m^ + 2q^ + iJmj^q'^ + m^q^ + mj^m'^ + g^ . (6.2) 

Then, using Eq. ( |A.1| ) and the results of Tab. ^ we can compute the threshold parameters. 
The first three of them are given in Tab. M for the four values of Scut investigated in our 
main fit. The uncertainties in Tab. ^ are obtained from a MC that takes into account all 
errors and correlations given in Tabs. ^ and ^. One should however note that the functional 
form of the threshold parameters, unlike A^ and A" , is determined by our model of S{s). 
Their values depend mainly upon the masses and widths of the resonances, most notably 
that of the K*{892). Since the pole of the K*{892) is very well determined in our fits, 
the uncertainties in the scattering lengths are accordingly small. Tab. ^ contains only the 
propagation of statistical uncertainties. The systematics uncertainty associated with the 
threshold parameters will be estimated in Sec. |^ 

7. Conclusions 

In this section we present our final results. They are obtained from the main fit displayed 
in Tab. y. Throughout this section, central values correspond to the average of the extrema 
found after the variation of Scut in Tab. y. Let us start with the mass and width of the 
K*{892) . To the statistical uncertainty one should add another source of error: the 
imperfect knowledge of the detector response. To that end, we rely on the original analysis 
performed by the Belle collaboration where it is found to be 0.44 MeV for the mass of the 
K*{892) and 1.0 MeV for its width^ [18]. In principle, one should include an uncertainty 
due to the residual dependence on Scut but Tab. ^ shows that the results are almost invariant 
under changes of this parameter. Therefore, this source can safely be neglected. Our final 
results for the mass and width of the i^*(892)^ defined from its pole position as in Eq. ( p.5| ) 
are then 

mK*i892)± = 892.03 ± (0.19)stat ± (0.44),ys MeV, 
ri^*(892)± = 46.53 ± (0.38)stat ± (l.O)sys MeV . (7.1) 



''The determination of the error due to detector effects is rather involved and depends on the model that 
is assumed for the analysis. Therefore, we take these values as mere estimates [50]. 
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Figure 3: Values for mass and width of the K*{892)^. On top, we show the PDG recommended 



values [41]. The other three results are obtained from the pole position as in Eq. (3^). The values 
for pole positions of Belle '07 [18] and Jamin et al. '08 [23] were computed in Ref. [24]. To the 
errors quoted in Ref. [24] we have added the systematics uncertainty discussed in the text. 



Before comparing this result with other analyses of the same data, one should note that 
in Refs. [18,23] a different definition of the mass of the K*[892) was used. Therefore, 
we have computed the pole position for the other analyses in order to harmonise the 
definition of mass. Moreover, the use of Eq. (|3.5D provides less model dependent results 
for the resonance parameters [36]. In Fig. |^, we compare the PDG recommended values 



m 



PDG 



ii-*(892)± = 891.66 ± 0.26 MeV and T^,fsg2)± = ^^-^ ^ ^-^ ^^'^ t^l] ^i^^ results for the 
mass and width of the K* (892)^ obtained from the pole position computed from the results 
of three different analyses of the Belle data set of r — )• Kirvr decays. Additional care should 
be taken when comparing these results since the PDG values are obtained chiefly from the 
parameters of Breit-Wigner-type expressions. On the basis of our results we claim that 
there is no tension between the mass found from r decays and the PDG recommended 
value provided the pole position prescription is used for the former. On the other hand, the 
PDG value for the width is only marginally compatible with the one from Eq. ([7.1|). The 
width from analyses of the Belle data on r — )■ KnUr tend to lower values. Let us conclude 
by quoting another unambiguous result that can be derived from our analysis: the point 
S7r/2 satisfying 6(3^^2) = '^/'^- Often, this point is used as the definition of the so-called 
visible or peak mass of a resonance since it is extracted from the direct comparison with 
experimental data^. In our fits, this value is also very stable with respect to changes in 
Scut and reads 

yi^T^ = 895.54 ± (0.01k,, MeV. (7.2) 



Our final values for A^ and A" come from the fits of Tab. ^ The results are again very 
stable with respect to changes in Scut- However, since now the statistical uncertainties are 
quite small, this model dependence contributes to the total error (specially in the case of 



®If the resonance in question is narrow and isolated enough from other resonances and, furthermore, if 
no background is present, then the mass and width obtained in this way should be equal to the pole mass 
definition [36]. We observe that the pole mass is close but not the same as the peak mass. This is due to 
the fact that we are not in the ideal situation stated before. 
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A'^). From the mean of values of Tab. y we obtain 

AV X 10^ = 25.49 ± (0.30)stat ± (0.06), 



A'l X 10^ = 12.22 ± (O.lO)stat ± (0.10),,,, . (7.3) 



Concerning A^, Fig. Q shows that the results from Ki^ and r — )• Kt^Vt decays are in very 
good agreement^. Our combined analysis produces a result in agreement with the others 
and with a rather small uncertainty. For A" the situation is somewhat different. Due to 
the restricted phase-space, quadratic fits of Ki^ data do not provide a good determination 
of A^ and dispersive analyses employ form factors with two subtractions, hence with only 
one subtraction constant determined directly from the data, namely A^. In Fig. |^, for 
Ki^ experiments, we display results derived from the two-times subtracted form factor of 
Ref . [26] . We compare our results also to an average of analyses that employ Eq. (1.1) for 



F+ [9]. Results from r decay data have a better precision, and are compatible with results 
from Ki^ experiments within their larger error bands. 

From the expansion of Eq (^]^) we can calculate the third coefficient of a Taylor series 
of the type of Eq. (^) as 

Scut 

X'l = a\ + ^a,a, + ml.^ f ds' ^ . (7.4) 

vr J [s'p 

Then, from the results of our fits, we find for A"' 

A'l' X 10^ = 8.87 ± (0.08)stat ± (0.05),,,, , (7.5) 

which is again compatible with the corresponding result of Eq. (5.7) of Ref. [24]. 

The ChPT expansion of F^{q^) at 0{p^) is governed by the low-energy constant Lg. 
Therefore, at this order, from our value of A| we can obtain Lg. It is not our aim here to 
carefully determine Lg , but it is certainly interesting to check the consistency of our results 
with the chiral expansion of F^. Using the 0{p'^) expressions of Ref. [58] with Fq = F^ we 
obtain 

Ll{mK*)\p2_p2 X 10^ = 5.19 ± (0.07)stat • (7.6) 

It is however well known that the dominant uncertainty is given by the truncation of the 
series at 0{p'^). As an estimate of 0{p^) effects we can employ Fq = F^^Fk which gives 

LlimK*)\p2_p p^ X 10^ = 6.29 ± (0.08)stat • (7.7) 

Our results agree with the one obtained in Ref. [59] from the pion electromagnetic form 
factor using ©(/) results: Ll{mK*) x 10^ = 5.70 ± 0.43. 

Our results for the phase space integrals of Ki^ decays have been collected in Tab. |5|. 
Following the procedure outlined above we find the final results given in the second column 



""For consistency we compare our results with dispersive analyses of Ki.^ decays. Data analyses that 



employ the quadratic Taylor expansion of Eq. (1.1) have much larger errors but agree as well with our 



numbers. For results from the quadratic form factor, see for instance [51] (ISTRA+), [52] (KLOE), [53] 
(NA48), and [54] (KTeV). 
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Figure 4: Results for A^ (left-hand panel) and A" (right-hand panel). The first four come from 
the dispersive representation of Ref. [26] fitted to Ki^ data of ISTRA+ [56], KTeV [27], KLOE [52], 
and NA48 [55]. A" of KLOE and NA48 are derived from their analyses in Ref. [57]. The results 
from Antonelli et al. '10. [9] are averages of fits using a quadratic form factor [see Eq. (1.1)]. The 
central portions of both panels display results from Moussallam '08 [28], Jamin et al. '08 [23], and 
our previous analysis [24]. 



of Tab. |7| In the same table, ■we display the results of the compendium performed in Ref. [9] 
from dispersive and quadratic fits to Ki^. Our results are compatible with those found in 
Ref. [9]. 





This Work 


Ki, disp. [9] 


Ki, quad. [9] 


^KO, 


0.15466(17) 


0.15476(18) 


0.15457(20) 


^n. 


0.10276(10) 


0.10253(16) 


0.10266(20) 


^K+ 


0.15903(17) 


0.15922(18) 


0.15894(21) 




0.10575(11) 


0.10559(17) 


0.10564(20) 



Table 7: Results for the Ki^^ phase-space integrals. Our results include a small uncertainty due to 
the dependence on s^ut- For comparison, we also give the results of Ref. [9] that come from averages 
of quadratic (quad.) and dispersive (disp.) analyses of X;, data. 



Finally, in Tab. ^ we present our final values for the i^vr P-wave 1=1/2 threshold 
parameters. These results are compared with other results found in the literature. Our final 
numbers include the statistical uncertainty as well as the (small) Scut dependence added 
in quadrature. Furthermore, we propagate the additional error of Eq. (7.1) in order to 
account for systematics. However, the precision obtained for the K*{S'd2) pole is such that 
our values have smaller uncertainties as compared to other determinations of the threshold 

1/2 

parameters. The main discrepancy observed is in the value of the effective range h^ that 
turns out substantially larger than that of Ref. [48] . In that reference, however, the authors 
already noted that their results could be affected by the uncertainties of LASS [46] data at 
energies above 1 GeV. The point where their phase equals it/2 is also shifted by 10 MeV as 
compared to ours. Therefore, since our data set is not contaminated with spurious strong 
interactions in the final state, we consider this discrepancy to be harmless. 

A final point concerning Kti interactions that should be address is the existence of 



16 



Table 8: Our final values for the threshold parameters compared with results found in the literature. 
In Rcf. [60] ChPT at ©(p'*) was used whereas in Ref. [61] ChPT at 0{p'^) was employed. Results 
from Ref. [62] are obtained within RChPT at O(p^) and in Ref. [48] a Roy-Steiner dispersive analysis 
of A'tt scattering was carried out. 

the controversial low-mass S'-wave isospin-1/2 resonance Kq(800) (or simply k,). In the 
description of the scalar form factor used here [31], a pole that can be identified with the 
K is present on the second Riemann sheet of the corresponding scattering amplitude [63] . 
Therefore, the success of our description of the spectrum in the low-energy region corrob- 
orates the existence of such a state. 

In conclusion, dispersion relations provide a technique to construct form factors valid 
for the description of r and kaon decay data. In the light of our results, we are confident 
that the use of dispersive form factors to fit the spectrum of r — )• KTrVr with restrictions 
from Ki^ experiments is a valid strategy towards the improvement of our knowledge of Kn 
form factors. Furthermore, some aspects of Kn dynamics can also be probed. New results 
for the spectrum of r — )• KttVt from other collaborations would offer a very good prospect 
to further improve our analysis. 
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A. Form factors 

A.l Vector form factor 

For the phase 6{s) needed in order to employ Eq. ( |2.3D we take a form inspired by the RChT 
treatment of Refs. [21,23] with two vector resonances. As described in greater detail in 
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Ref. [24], 6{s) can be cast into the following form 

Im /4 



6{s) = tan 



Re /+(s) 



(A.l) 






where 



The first term in the right-hand side of Eq. ( |A.2| ) corresponds to the K*(892) whereas the 
second represents the contribution of the second vector resonance i^*(1410). The mixing 
parameter 7 is obtained from the fits and Hktt{s) is the one-loop Kit bubble integral, 
whose precise definition is given in Refs. [21,58]. The denominators D{mK*,'yK*) are 

-D(m„,7„) = ml - s - KnRe Hk-k{s) -irhnlnis) , (A.3) 

where the constants 

„„ . 1^^^^ (A.4) 

are defined so that — z«;„Im Hk-k{s) = —imn^nis) and the running width of a vector 
resonance is taken to be 

The phase-space function aKwis) is given by crKn{s) = 2 qK-K{s)/y/s, whereas qk-k is defined 
in Eq. ( |3.2D . The model parameters m„ and 7„ are not the physical resonance mass and 
width. Physical values are obtained solving the equation D{rhn-,'7n) = for complex values 
of s. Consequently, our definition of physical mass and width is given by Eq. ( |3.5| ). 

A. 2 Scalar form factor 

The procedure adopted in Ref. [29] is to solve the multi-channel Muskelishivili-Omnes 
problem for 3 channels (where 1 = KTr,2 = Krj and 3 = Kt]'). Each of the scalar form 
factors Fq , where k represents the channel, is then coupled to the others via 



FoHs) 



It ["'"-m^S^^ (A.«) 



vr ^-^ I (s' — s — ie) 



In the last equation, Sj is the threshold for channel j, (Tj{s) are two-body phase-space 
factors and tg '' are partial wave T-matrix elements for the scattering k ^- j. The form 
factors are obtained solving the coupled dispersion relations arising from Eq. ( |A.6| ). This is 
done imposing chiral symmetry constraints and using T-matrix elements from Ref. [63] that 
provide a good description of scattering data. Within the elastic approximation, Eq. ([A. 61 ) 
reduces to the usual single-channel Omnes equation. 
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